Synchronization of Multiple Coupled rf-SQUID Flux Qubits 
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A practical strategy for synchronizing the properties of compound Josephson junction rf-SQUID 
qubits on a multiqubit chip has been demonstrated. The impacts of small (~ 1%) fabrication varia- 
tions in qubit inductance and critical current can be minimized by the application of a custom tuned 
flux offset to the CJJ structure of each qubit. This strategy allows for simultaneous synchronization 
of the qubit persistent current and tunnel splitting over a range of external bias parameters that is 
relevant for the implementation of an adiabatic quantum processor. 
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Despite daunting obstacles, there is considerable in- 
terest in the development of solid state quantum infor- 
mation processors. This interest is fueled by the hope 
that breakthroughs in device fabrication will eventually 
facilitate the realization of large scale quantum proces- 
sors whose performance could surpass that of classical 
computers. Implementations based upon superconduct- 
ing qubits have received particular attention [1]. Con- 
siderable effort has been made in studying noise in such 
circuits [2, 3]. An equally pressing matter is fabrication 
variability as qubits are acutely sensitive to variations 
in device parameters [4, 5]. Current state of the art su- 
perconducting fabrication technology is limited, at best, 
to ~ 1% spreads in parameters such as Josephson junc- 
tion critical currents and qubit inductances. The extent 
to which this variability affects the performance of a su- 
perconducting quantum computer is an important open 
problem. Therefore, it is relevant to demonstrate super- 
conducting qubit designs and methods of operation that 
are insensitive to small variations in device parameters. 

We wish to focus on a quantum Ising spin glass simu- 
lator [6, 7] constructed from superconducting flux qubits 
[8] . Such a device could be useful for solving optimization 
problems [9]. Each qubit i serves as a spin-1/2 subjected 
to transverse and longitudinal biases A, and e-s = [liBi, 
respectively. Here, \ii represents the effective magnetic 
moment and Bi an externally controlled magnetic held. 
Pairwise couplings are realized by Jjj = MijUifij, where 
Mi_j is an externally controlled parameter. The system 
Hamiltonian at any time during operation has the form 
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A particular adiabatic quantum algorithm, such as that 
described in Ref. [6], may require all fa and Aj to be 
nominally equivalent between qubits. While this choice 
of algorithm is by no means unique, it does represent 
the simplest implementation of an optimization proce- 




FIG. 1: (Color online) Schematic diagram of a CJJ rf-SQUID. 



dure that utilizes quantum adiabatic evolution. The ob- 
jective of the work presented herein was to develop a 
practical strategy for minimizing the differences in qubit 
parameters between superconducting flux qubits due to 
fabrication variations via in-situ tunable device biases. 

One useful implementation of a superconducting flux 
qubit is the compound Josephson junction (CJJ) rf- 
SQUID [5], as depicted in Fig. 1. Here, a main loop of 
superconducting wire of inductance L q is interrupted by a 
smaller loop of inductance L c jj with two Josephson junc- 
tions of critical current I c \ and I C 2 ■ The CJJ and main 
loop are subjected to external fluxes = ^ofcj}/^ n 
and $ q = $> (p*/2n, respectively ($ = h/2e). The 
Hamiltonian for this system can be written as 
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where Cj and Eji — J C i<I>o/27r represent the capaci- 
tance and Josephson energy of junction i, respectively, 
and [$oVi/27r, Qj] = ihSij. The inductive terms origi- 
nate from the two closed loops with n G {q, cjj}, L q = 
L + £ C jj/4 and U n = ($ /27r) 2 /L„. The qubit and 
CJJ loop phases are defined as tp q = (tp\ + if2) /2 and 
ifcjj = ipi — (f2, respectively. This 2-dimensional system 
can be reduced to an effective 1-dimensional Hamiltonian 
if L q S> L c]i because the plasma energy of the CJJ loop 
will be much higher than that of the main rf-SQUID loop. 
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Setting ifcjj = ip%-! and combining the Josephson terms, 
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where C p = C\ + C2, [^of q /2Tr, Q q ) = i/l and (3± = 
27rL g (7 c i ± I C 2) /$o- Focussing upon the two lowest ly- 
ing states in the regime /3 e fi < — 1 (7r phase shifted 
flux qubit), one can recast Eq. 3 as a qubit Hamilto- 
nian H q = -\ [e<7 z + Aa x ], where e = 2 - *°). 

Denoting the ground and first excited state of Eq. 3 
at 4>;£ = $° by |+) and |— ), respectively, the spin 
states can be expressed as ||) = (|+) + |— )) /y/2 and 
I J,) = (|+) — |— )) /V2- The persistent current is then 
defined by |7|| = |(|| ($ g - $°) /L g ||>|. The tunneling 
energy is given by A = (— | H |— ) - (+| 7Y |+). 

The CJJ rf-SQUID was first proposed as a means of 
providing in-situ tuning of the tunnel barrier U q f3 e ft via 
[5]. By allowing for small relative offsets in the 
CJJ bias <5cI> C jj , it is possible to simultaneously minimize 
differences in \l q \ and A between qubits with slightly 
different L q and I c = I c \ + I C 2, thus synchronizing 
their properties. Consider |/^| and A in the regime 
UqPefi S> h\o p = h/y/LqCp. In this scenario, is 
primarily determined by the position of the minima of 
V(ifq) with only an extremely weak dependence upon 
C p . In order to maintain constant \l q \ in the pres- 
ence of small variations in a g {L q ,I c \ 1 the condition 
is /3 e ff(a, $cij) = A=ff(a + $a,®*y } + <5$ C jj)- To first or- 
der in 6a/ a, <5<f> cj j ps (^o/ 71 ") [ cot (^cjj/^o)] <W a - For 
|5a/a| = 0.05, one obtains |5$ c jj| ~ 15m$ - Further- 
more, one can use the WKB approximation [10] to write 
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where ±a represent the classical turning points straddling 
the local maximum in V(ip q ). The resultant form for A 
reveals that A(a, = A (a + <5a, + <5<f> c jj), where 
<5$ cjj ps 7($ /tt) [cot (7r<&^j/<E>o)] 5a/a, with 7 ~ 1 for 
a = L q , I c and 7 <C 1 for a — C p . Interestingly, A shows 
a relatively weak dependence upon C p as compared to L q 
and I c . Thus perturbations of L q and I c (< 5%) result in 
approximately the same shift of the CJJ bias-dependence 
of both I J* I and A. In contrast, perturbations of C p 
(< 5%) have negligible impact upon \l q \ but do influence 
the CJJ bias dependence of A. 
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FIG. 2: (Color online) Schematic of a portion of a multiqubit 
chip, bias line configuration and room temperature electron- 
ics. Ellipses indicate devices (alternating between qubit and 
coupler) that extend beyond the scope of the diagram. 



The above observations indicate that one can compen- 
sate for small variations of L q and I c between CJJ rf- 
SQUID qubits by the application of custom tuned CJJ 
bias offsets. For typical device parameters, A/h varies 
from - 1MHz to - 10 GHz for 1 < |/3 off | < 1.3: If the 
qubit has been designed with j3 + > 1.5, then the range of 
that is relevant for operation will be < 50 m$o wide. 
Thus, one can choose a unique reference CJJ bias <E>°jj 



in the center of the operating regime for each qubit such 
that |/|| — and A(<I>^ jj - are synchronized. 

In order to experimentally assess the CJJ synchroniza- 
tion strategy we focussed on a subset of CJJ rf-SQUID 
flux qubits embedded in a larger lattice of such devices 
(see Fig. 2). Each qubit is connected to three others 
via in-situ tunable monostable CJJ rf-SQUID couplers, 
which we treat as classical mutual inductances [11]. We 
isolated a linear chain of six qubits by setting the inter- 
vening couplers to maximum antifcrromagnctic coupling 
and the remaining unused couplers to provide zero cou- 
pling. Each qubit 's state was probed via a dedicated 
dc-SQUID magnetometer. The chip was fabricated from 
an oxidized Si wafer with Nb/Al/AlaOa/Nb trilayer junc- 
tions and three Nb wiring layers separated by sputtered 
Si02- It was mounted to the mixing chamber of a dilution 
refrigerator regulated at T = 35 mK inside a PbSn super- 
conducting magnetic shield with a residual field in the 
vicinity of the chip < 9 nT. External current biases were 
provided by a custom-built programmable room temper- 
ature 128-channel current DAC. Low pass filters (LPFs) 
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FIG. 3: (Color online) CJJ and flux bias waveforms versus 
time t for source (solid) and detector (dashed) qubits. 



with f c ~5 MHz were constructed from a combination of 
lumped element and copper powder filters secured to the 
mixing chamber. All mutual inductances and residual 
flux offsets were calibrated in-situ. 

To measure \lq \ we first employed dc-SQUIDs as mag- 
netometers. To begin, we initialized a given qubit, here- 
after referred to as the source qubit, in the state |T), 
raised U q j3 e g to maximum height (<£>^j = —$0) and mea- 
sured the change in flux sensed by its dc-SQUID. The 
process was repeated for the qubit initialized in the state 
IJ.) and the difference between the two measurements 
recorded. Knowing the readout to qubit mutual induc- 
tance M ro _ q = 6.46±0.17pH, obtained from an indepen- 



Pcrforming the measurement for both signs of <& s q and 
taking the difference between the two resultant values of 
<f>^ yielded 2M c g Given J£ we then inferred 

1 J I 9 1 I 1 I max 

M e g = 1.35 ± 0.04 pH for the 5 intervening couplers in 
the chain of 6 qubits. ft was then possible to scale maps 
of 2M eS \IP\ versus ^ to extract 

To measure A of each qubit we used two methods: In 
the incoherent regime one can utilize macroscopic res- 
onant tunneling (MRT) to trace out decay rate curves 
and extract A from fitting parameters [12]. The range 
of A that could be probed by this method had a practi- 
cal upper bound on account of the relatively low band- 
width of our bias lines. In the coherent regime, we em- 
ployed a 2-qubit procedure involving the waveform pat- 
tern shown in Fig. 3 in which $^ was scanned through 
the domain [—3, +3] m<& and <fr d was again adjusted 
via a software feedback procedure to determine the shift 
in detector qubit degeneracy point at each & q . Know- 
ing |/p| ($cjj) allowed one to convert and & d into 
ei and e 2 , respectively. For two coupled qubits in the 
limit A2 — > the eigenenergies of Eq. 1 are given by 

' e a l, 



E 1± = I[±F(-l)-e 2 ] 
where F(x) = \J (ei 



and E 2± = i [±F(+1) ■ 

Using Boltzmann 
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statistics, one can calculate the particular bias e 2 = e 2 
for which the detector qubit will be found with Pf = 1/2: 

. F( + l)-F(-l) / l + e -F(+l)/*BT N 

= + fc B Tln [ - : _ F ,_ 1W ^„ T , ) (5) 
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dent measurement, we determined \If,\ = \If, = 

_ _ I H I max I h I v C JJ 

— $0) for each qubit. This measurement technique pro- 



Eq. 5 reduces to 
\g), with \g) rep- 



IP 
9 I 



vided reliable results only if Af ro _ q 
£$ ro ~ 2 m$o represents a flux resolution limit imposed 
by the width of the dc-SQUID switching current distri- 
bution. In order to clearly resolve |/p| with U q /3 e g sup- 
pressed, we utilized a second qubit, hereafter referred to 
the detector qubit, that was coupled to the source qubit 
via a coupler with effective mutual inductance M e g. Re- 
ferring to Fig. 3, the sequence began with both qubit 
V(ip q ) monostable = — $o/2) and biased to their 

degeneracy points (Q q — $° = 0) (i). Next, the source 
qubit was partially annealed to an intermediate CJJ bias 
— $0 < &cjj < —^0/2 in the presence of a small bias 
§ s q = ±2.1 m$o in order to initialize its state (ii). There- 
after, the detector qubit was fully annealed ramped 
to — $0) in the presence of a variable bias Q d (iii). Fi- 
nally, the source qubit was fully annealed (iv) , both qubit 
flux biases are returned to their degeneracy point (v) and 
the state of the detector qubit is read (not shown). This 
annealing cycle was embedded inside a software feedback 
loop which adjusted <& d until the particular bias for which 
the detector qubit could be found in ||) with probabil- 
ity Pj = 1/2 was found to within a specified precision. 



> 5<& TO , where T, one can fit traces of e% 



resenting the groundstate of the source qubit. Given in- 
dependent calibrations of for both qubits, M e g and 

versus ei with Eq. 5 to ex- 
tract Ai. This procedure is similar in spirit to that of 
Ref. [13]. In practice, the 2-qubit method was found to 
be reliable only if A > T and M ef r|/P| > <5<f>„, where 
<5<f> n represents an rms low frequency flux noise experi- 
enced by the detector qubit. These constraints imposed 
lower and upper bounds, respectively, upon the range of 
Ai that could be probed via this latter method. 

Measurements of the CJJ bias dependence of and 
A are shown in Fig. 4. Here, we have shifted the CJJ bias 
for each qubit by a unique <f>®-^ (as summarized in Table I) 
chosen such that A($°jj)//i = 10 7 Hz. It can be seen that 
the six sets of |/p| data lie atop one another to within 
the measurement uncertainty over the range of — <3>|?jj 
for which A varies by five orders of magnitude. The A 
data show reasonable synchronization, albeit the results 
for q 2 show higher A in the coherent regime and slightly 
faster exponential decay as a function of at small A. 
Otherwise, the values of A from the other 5 qubits are 
synchronized to within 20% over the range of — "J 1 ^ 
shown. The 2-qubit method for extracting A proved par- 
ticularly susceptible to corruption by low frequency flux 
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FIG. 4: (Color online) a) and b) A//i as a function of 
synchronized CJJ bias. A/h from the 2-qubit and MRT mea- 
surement procedure are denoted as (g \ Ig \g) and MRT in the 
legend, respectively. Solid curves are theoretical predictions 
using the mean device parameters quoted in Table I. 



noise. For these qubits, drift measurements of the type 
reported in Ref. [3] revealed 1// noise spectral densities 
with a mean amplitude y/S$(l Hz) = 14 ± 2 / v44z. 
Efforts to refine the 2-qubit method and to reduce 1/f 
noise in our devices are ongoing. 

The results of simultaneously fitting an< ^ 
A($^jj) using the eigenstates of Eq. 3 for each qubit are 
summarized in Table I. The quality of the fits proved 
most sensitive to I c and comparably less sensitive to the 
choice of L q and C p . It is probable that fabrication vari- 
ations between Josephson junctions, roughly ±1% of the 
target I Cl are the prime source of inter-qubit variability 
on this particular chip. The theoretical predictions for 
rq \ (^cjj) an d ^(^cjj) usm g the mean device parameters 
listed at the bottom of Table I are shown in Fig. 4. 



Qubit 


$cjj (m$o) 


Lq (PH) 


C p (fF) 


ic (M) 


1 


-789 ± 5 


200 ±2 


56 ±1 


2.58 ±0.01 


2 


-774 ± 5 


202 ±2 


56 ±1 


2.65 ±0.01 


3 


-781 ± 5 


200 ±2 


57 ±1 


2.63 ±0.01 


4 


-784 ± 5 


202 ±2 


55 ±1 


2.59 ±0.01 


5 


-777 ± 5 


200 ±2 


56 ±1 


2.65 ±0.01 


6 


-785 ± 5 


202 ±2 


54 ±1 


2.59 ±0.01 


Mean 


-782 ± 12 


201 ± 1 


56 ±1 


2.62 ±0.03 



TABLE I: Relative CJJ bias shifts 4?^ and device parameters 
obtained by simultaneously fitting ($cjj) an d ^(^cjj)- 



Conclusions: A method for synchronizing the proper- 
ties of multiple coupled CJJ rf-SQUID flux qubits with 
a small spread of device parameters due to fabrication 
variations has been demonstrated. Both theory and ex- 
periment indicate that the application of a custom-tuned 
flux bias to each qubit CJJ loop is sufficient to compen- 
sate for ±1% differences in critical current. This strategy 
may prove to be an important step in the development 
of practical adiabatic quantum information processors. 

We thank J. Hilton, G. Rose, P. Spear, A. Tcaciuc, 
F. Cioata, E. Chappie, C. Rich, C. Enderud, B. Wilson, 
M. Thorn, S. Uchaikin and M.H.S. Amin. Samples were 
fabricated by the Microelectronics Laboratory of the Jet 
Propulsion Laboratory, operated by the California Insti- 
tute of Technology under a contract with NASA. S.Han 
was supported in part by NSF Grant No. DMR-0325551. 



[5] 



[»] 



Electronic address: rharris@dwavesys.com 
URL: www.dwavesys.com 

John Clarke and Frank K. Wilhelm, Nature 453, 1031 
(2008), and references therein. 

R.W. Simmonds et a/., Phys. Rev. Lett. 93, 077003 

(2004) ; J.M. Martinis et ai, Phys. Rev. Lett. 95, 210503 

(2005) ; R.C. Bialczak et ai, Phys. Rev. Lett. 99, 187006 

(2007) ; O. Astafiev, Y.A. Pashkin,Y. Nakamura, T. Ya- 
mamoto and J.S. Tsai, Phys. Rev. Lett. 93, 267007 
(2004); J. A. Schreier et ai, Phys. Rev. B 77, 180502(R) 

(2008) . 

T. Lanting et ai, Phys. Rev. B 79, 060509(R) (2009). 

R.C. Ramos et ai, IEEE Trans. Appl. Supercond. 11, 

998 (2001); Jens Koch et ai, Phys. Rev. A 76, 042319 

(2007); J.E. Mooij et ai, Science 13, 285 (1999); F.G. 

Paauw, A. Fedorov, C.J. P.M. Harmans and J.E. Mooij, 

Phys. Rev. Lett. 102, 090501 (2009). 

S. Han, J. Lapointe and J.E. Lukens, Phys. Rev. Lett. 

63, 1712 (1989); S. Han, J. Lapointe and J.E. Lukens, 

Phys. Rev. Lett. 66, 810 (1991). 

E. Farhi et ai, Science 292, 472 (2001); 

D. Aharonov, W. van Dam, J. Kempe, Z. Landau, and 

S. Lloyd, SIAM Journal of Computing 37,166 (2007); 

Jacob D. Biamonte and Peter J. Love, Phys. Rev. A 78, 

012352 (2008); Jingfu Zhang et ai, Phys. Rev. A 79, 

012305 (2009). 

W.M. Kaminsky and S. Lloyd, in Quantum Comput- 



5 



ing and Quantum Bits in Mesoscopic Systems, MQC 2 
(Kluwer Academic, New York USA, 2003). 
[9] E. Boros, P.L. Hammer and G. Tavares, J. Heuristics 13, 
99 (2007). 

[10] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 

Vol. 3, 3rd ed., 179, (Butterworth-Heinemann, 1977). 
[11] A. Maassen van den Brink, A.J. Berkley, and 



M. Yalowsky, New J. Phys. 7, 230 (2005). 
[12] M.H.S. Amin and D.V. Averin, Phys. Rev. Lett. 100, 
197001 (2008); R. Harris et al, Phys. Rev. Lett. 101, 
117003 (2008), 

[13] A. Izmalkov et al, Phys. Rev. Lett. 101, 017003 (2008). 



